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Abstract. The decay rates and spectroscopy of the QQ (Q € c, b) mesons are 
computed in non-relativistic phenomenological quark antiquark potential of the type 
V(r) — + Ar v , (CPP„) with different choices v. Numerical solution of the 

schrodinger equation has been used to obtain the spectroscopy of QQ mesons. The spin 
hyperfine, spin-orbit and tensor components of the one gluon exchange interaction are 
employed to compute the spectroscopy of the few lower S and orbital excited states. 
The numerically obtained radial solutions are employed to obtain the decay constant, 
di-gamma and di-leptonic decay widths. The decay widths are determined with and 
without radiative corrections. Present results are compared with other potential model 
predictions as well as with the known experimental values. 
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1. Introduction 

Heavy flavour hadrons play an important role in several high energy experiments as 
well as in the understanding of the theories like QCD, NRQCD, pNRQCD, vNRQCD 
and effective field theories. The BES at the Beijing Electron Positron Collider (BEPC), 
E835 at Fermilab, and CLEO at the Cornell Electron Storage Ring (CESR) are able 
to collect the huge data on heavy flavour mesons. Where as B-meson factories, BaBar 
at PEP-II and Belle at KEKB are working on the observation of new and possibly ex- 
otics quarkonia states. The CDF and experiments at Fermilab measuring heavy 
quarkonia production from gluon-gluon fusion in pp annihilations at 2 TeV^. Also some 
other experiments like ZEUS and HI at DESY are studying charmonia production in 
photon-gluon fusion. The study related to the charmonia production and suppression 
in heavy- ion collisions are being looked by PHENIX, STAR and NA60. All these exper- 
iments are capable of observing new states, new production mechanisms, new decays 
and transitions, and in general to the collection of high statistics and precision data 
sample. In the near future, even larger data samples are expected from the BES-III 
upgraded experiment, while the B factories and the Fermilab Tevatron will continue to 
supply valuable data for few years. Later on, the LHC experiments at CERN, Panda at 
GSI etc are capable of offering future opportunities and challenges in this field of heavy 
flavour physics [lj. 

On the theoretical side, heavy quarkonium provides testing and the validity of pertur- 
bative QCD, potential models and lattice QCD calculations [2|. The investigation of 
the properties of mesons composed of a heavy quark and antiquark (cc, be, hh) gives 
very important insight into heavy quark dynamics and to the understanding of the con- 
stituent quark masses. The theoretical predictions of the heavy quarkonia cc, he and 
bb mesons have rich spectroscopy with many narrow states of charmonium lying un- 
der the threshold of open charm production (3j Hj and of botomonium lying under the 
threshold of B — B production. Many of these states have not confirmed or understood 
by experiments [5]. However, there have been renewed interest in the spectroscopy of 
the heavy flavoured hadrons due to number of experimental facilities (CLEO, DELPHI, 
Belle, BaBar, LHCb etc) which have been continuously providing and expected to pro- 
vide more accurate and new informations about these states at the heavy flavour sector. 
At the hadronic scale the non-perturbative effects connected with complicated struc- 
ture of QCD vacuum necessarily play an important role. All this leads to a theoretical 
uncertainty in the QQ potential at large and intermediate distances. It is just in this 
region of large and intermediate distance that most of the basic hadron resonances are 
formed. So the success of theoretical model predictions of most of the hadronic prop- 
erties with experiments can provide important information about the quark-antiquark 
interactions. Such information is of great interest, as it is not possible to obtain the QQ 
potential starting from the basic principle of the quantum chromodynamics (QCD) at 
the hadronic scale. 

Among many theoretical attempts or approaches to explain the hadron properties based 
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on its quark structure very few were successful in predicting the hadronic properties 
starting from its spectroscopy to decay rates. The nonrelativistic potential models with 
Buchmiiller and Tye [6], Martin [TJ [8J E], Log [10, E], Cornell [12] etc. were successful in 
predictions of the spectra of the heavy flavour mesons while the Bethe-Salpeter approach 
under harmonic confinement [13] were successful at low flavour sector. Though there 
exist relativistic approaches for the study of the different hadronic properties [TH[T5| [T6]. 
the non-relativistic models have also been equally successful at the heavy flavour sector. 
For the theoretical predictions of different decay rates most of the models require sup- 
plementary corrections such as higher order QCD effects, radiative contributions etc. 
Even in some cases rescaling of the model radial wave functions are also being con- 
sidered. However the NRQCD formalism provides a systematic approach to study the 
decay properties like the di-gamma and the di-lepton decays. These partial decay widths 
provide an account of the compactness of the qurkonium system which is an useful in- 
formation complementary to spectroscopy [17] . Thus, in this paper we make an attempt 
to study the properties like mass spectrum, decay constants and other decay properties 
of the QQ systems (Q e b, c) based on a phenomenological coulomb plus power potential 
(CPP V ). Here, we consider different choices of the potential power index v to study the 
properties of the mesonic systems upto few excited states. 

2. Nonrelativistic Treatment for Heavy Quarks 

There are many theoretical approaches both relativistic and nonrelativistic to study 
the heavy quark systems [niDBDIlIIlEIlEaE^ However 
their predictions suggests the successes of the nonrelativistic treatment for the heavy 
flavour quark-antiquark system [TT]. The relativistic invariant theory for example light- 
front QCD [29] though deals with different aspects of QCD, under non-relativistic 
approximations, reproduces the results comparable to the non-relativistic quark- 
potential models [29]. In the center of mass frame of the heavy quark-antiquark 
system, the momenta of the quark and antiquark are dominated by their rest mass 
m Q,Q ^ ^qcd ~ \p |) which constitutes the basis of the non-relativistic treatment. For 
examples NRQCD formalism for the heavy quarkonia, the velocity of heavy quark is 
chosen as the expansion parameter |30j . 

Hence, for the study of heavy-heavy bound state systems such as cc, be and bb, we 
consider a nonrelativistic Hamiltonian given by [2H (25J [26] 



momentum of each quark and V(r) is the quark antiquark potential. Though linear plus 
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Table 1. The model Parameters employed in the present study: (The potential 
strength A for different power index, v is given in GeV !y+1 ) 



V 


A (cc) 


A (6c) 


A (66) 


0.5 


0.3630 


0.4085 


0.4600 


0.7 


0.3034 


0.3582 


0.4430 


0.8 


0.2784 


0.3366 


0.4358 


0.9 


0.2559 


0.3169 


0.4296 


1.0 


0.2355 


0.2986 


0.4237 


1.1 


0.2170 


0.2817 


0.4180 


1.3 


0.1846 


0.2513 


0.4080 


1.5 


0.1573 


0.2246 


0.3984 



a c (cc) = 0.40, a c {bc) = 0.34, a c {bb) = 0.30, 
m c = 1.24 GeV and m h = 4.50 GeV 



coulomb potential is a successful well studied non-relativistic model for heavy flavour 
sector, their predictions for decay widths are not satisfactory owing to the improper 
value of the radial wave function at the origin compared to other models [25]. Thus, in 
the present study we consider a general power potential with color coulomb term of the 
form 

V(r) = -^ + Ar» (3) 
r 

as the static quark-antiquark interaction potential. This potential belong to the special 

choices of the generality of the potentials, V(r) = —Cr a + Dr 13 + Vq [3H E21 [33] with 

V = a = —1, (3 = u. For the present study, the power index range of 0.1 < v < 2.0 

4 

have been explored. Here, for mesons, a c = -a s , a s being the strong running coupling 

3 

constant, A is the potential parameter similar to the string strength and v is a general 
power, such that the choice, v — 1 corresponds to the coulomb plus linear potential. 
The different choices of v here, correspond to different potential forms. In general, the 
potential parameter A can also be different numerically and dimensionally for each 
choices of v. In the present study of heavy-heavy flavour mesons, we employ the 
numerical approach |34J to generate the Schrodinger mass spectra. 



3. Spin-Dependent forces in QQ States 

In general, the quark-antiquark bound states are represented by n 2S+1 Lj, identified 
with the J PC values, with J = L + S, S = Sq + Sq, parity P = (— 1) L+1 and the 
charge conjugation C = (— 1) L+S with (n, L) being the radial quantum numbers. So 
the S-w&ve (L = 0) bound states are represented by J PC = 0~ + and 1 respectively. 
The P-wave (L = 1) states are represented by J PC = l + ~ with L = 1 and S = while 
J PC = ++ , 1 ++ and 2 ++ correspond to L — 1 and 5=1 respectively. Accordingly, 
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Figure 1. Behavior of A with the potential index v for different QQ systems 



the spin-spin interaction among the constituent quarks provides the mass splitting of 
J = 0~ + and 1 states, while the spin-orbit interaction provides the mass splitting of 
J PC = ++ , 1 ++ and 2 ++ states. The J = l +_ state with L = 1 and S = represents 
the center of weight mass of the P-state as its spin-orbit contribution becomes zero, while 
the two J = l + ~ singlet and the J = 1 ++ of the triplet P-states form a mixed state. 
The Z)-wave (L = 2) states are represented by J PC = 2~ + with L = 1 and S = while 
J PC = 3 , 2~ + and 1 h correspond to L = 2 and S = 1 respectively. 
For computing the mass difference between these states, we consider the spin dependent 
part of the usual OGEP given by [35] as 



and 

a c (3(5 ? -n)(5 ? -n) - S- S) _ f 

VT\T) = TT o 1 71= - 5) 

M Q m Q r 3 ' r 

The spin average mass for the ground state is computed for the different choices of v 
in the range, 0.5 < v < 1.5. The model parameters used here are listed in Table [U 
The potential parameter A are fixed for each choices of v so as to get the experimental 
ground state spin average masses of QQ systems. The spin average masses of cc is 
computed using the experimental ground state mass of M Vc = 2.980 GeV and Mj/^ = 
3.097 GeV [5], while the experimental values of Mx = 9.460 GeV and theoretically 
predicted values for rjb, M Vb = 9.400 GeV [15] are used to get the centre of weight mass 
of bb system. For the be meson we use the experimental mass of Mb c = 6.286 GeV [5] 
and the theoretically predicted value of Mb* = 6.332 GeV pj[|. For the nJ state, we 
compute the spin-average or the center of weight mass from the respective experimental 
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values as 




Ej2(2J+l) M nJ 
Ej2(2J+l) 



(6) 



In the case of quarkonia (cc and bb systems) many orbital excited states are known. 
Theoretical predictions of all these states and their decay widths are also being studied. 
But in many cases, the decay widths and the spin splitting between different J values 
are not well reproduced. Both the decay widths and the level splitting of the spectra 
due to the one gluon exchange interaction terms are related to the values of the radial 
wave function or its derivatives at the origin. Thus, the inappropriate description of the 
QQ radial wave function led to the disparity among the different model predictions of 
the decay widths and level splitting. In some cases, for better predictions of the excited 
spectra of quarkonia, the strong running coupling constant a s are evaluated interms of 
the average kinetic energy of the quark-antiquark pair at a given state. Accordingly, 
different excited states corresponds to have different values of a s [35], [36] . However, the 
radial wave functions are found to be less sensitive to the changes in a s compared to 
similar changes in the values of the strength of the confining part of the potential. Hence, 
in this paper, we allow A to vary mildly with radial quantum number (n = 0, 1, 2...) as 
A = — — — j- . The variation in A can be justified by similar arguments for the changes in 



as with the average kinetic energy. Here, as the system get excited, the average kinetic 
energy increases and hence the potential strength (the spring tension) reduces. With 
this mild state dependence on the potential parameter A, we obtain the excited spectra 
as well as the right behavior for the radial wave functions . The computed values of the 
radial wave function at the origin for (n + 1)S states are listed in Table [2] for all the QQ 
combinations. Using the spin dependent potential given by Eqn. H] and El we compute 
the masses of the different n 2S+1 Lj states of cc, be and bb mesons. Better stastics with 
respect to the experimental values are observed with our predictions of these states for 
the potential index lying between 0.7 to 1.3. Thus we list our predicted properties in 
this range of potential index only. The computed masses of the QQ mesonic states are 
listed in Table [3] in the case of cc , in Table H] in the case of be and in Table [5] in the 
case of bb systems along with the available experimental values as well as other model 
predictions. Fig [1] shows the behavior of A with the potential index v that provide us 
the ground state center of weight masses for all the three (cc, be and bb) combinations 
of QQ systems. 

4. The Decay constants of the heavy flavoured mesons 

The decay constants of mesons are important parameters in the study of leptonic or non- 
leptonic weak decay processes. The decay constants of pseudoscalar (f P ) and vector (f v ) 
states are obtained by prarameterizing the matrix elements of weak current between the 
corresponding mesons and the vacuum as 



(n+l)2 



(0\Qr l5 Q\P,(k)) = if P k» 



(7) 
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Table 2. The radial Wave function \R ns (0)\ 2 (in Get/ 3 ) of QQ systems m various 
potential models including CPP„. 



Mesonic 


Potential 




IS 


2S 




3S 


AS 


55 


6S 


System 


TV T 11 

Model 
















cc 


(GPF v ),v = 


= 0.5 


0.420 


0.198 


0.136 


0.106 


0.088 


0.075 






0.7 


0.529 


0.295 


0.221 


0.182 


0.158 


0.141 






0.8 


0.577 


0.347 


0.270 


0.229 


0.202 


0.183 






0.9 


0.622 


0.400 


0.323 


0.280 


0.252 


0.232 






1.0 


O.oo2 


0.454 


0.379 


0.337 


0.309 


0.288 






1.1 


0.700 


0.509 


0.439 


0.399 


0.372 


0.352 






1.3 


0.767 


0.623 


0.569 


0.538 


0.517 


0.502 




TV T J_ * II-7I 

Martin [7] 




0.979 


0.545 


0.390 


0.309 


0.257 


0.222 




Log [TD] 




0.796 


0.406 


0.277 


0.211 


0.172 


0.145 




Cornell [37] 




1.458 


0.930 


0.793 


0.725 


0.683 


0.654 




T— J 1 11 

Bucnmuller- 


-Tye [B] 


0.794 


0.517 


0.441 


0.404 


0.381 


0.365 




T • 1 1 1 

Licntenberg 


-Wills [38] 


1.121 


0.693 


0.563 


0.496 


0.453 


0.423 


be 


(GPP v ),u = 


- 0.5 


0.886 


0.411 


0.280 


0.217 


6.865 


0.154 






0.7 


1.109 


0.609 


0.454 


0.373 


7.221 


0.288 






0.8 


1.207 


0.714 


0.553 


0.468 


7.414 


0.374 






0.9 


1.298 


0.823 


0.661 


0.573 


7.615 


0.473 






1.0 


1.381 


0.933 


0.776 


0.688 


7.823 


0.587 






1.1 


1.457 


1.047 


0.898 


0.814 


8.039 


0.716 






1.3 


1.594 


1.278 


1.163 


1.098 


8.487 


1.020 




Martin [7] 




1.720 


0.957 


0.685 


0.452 


0.452 


0.390 




Log pig 




1.508 


0.770 


0.524 


0.401 


0.325 


0.275 




s~H ii 1 T 1 

Cornell [37] 




3.191 


1.769 


1.449 


1.297 


1.205 


1.141 




T ~\ 1 11 

Bucnmuller- 


-Tye [6] 


1.603 


0.953 


0.785 


0.705 


0.658 


0.625 




Lichtenberg 


-Wills [38] 


2.128 


1.231 


0.975 


0.846 


0.766 


0.711 


bb 


(CPP„),z/ = 


= 0.5 


4.222 


1.750 


1.151 


0.876 


0.716 


0.610 






0.7 


5.101 


2.534 


1.828 


1.479 


1.265 


1.118 






0.8 


5.487 


2.948 


2.214 


1.840 


1.607 


1.443 






0.9 


5.843 


3.377 


2.632 


2.244 


1.996 


1.820 






1.0 


6.170 


3.814 


3.078 


2.687 


2.433 


2.251 






1.1 


6.470 


4.259 


3.551 


3.168 


2.917 


2.735 






1.3 


7.006 


5.173 


4.575 


4.249 


4.034 


3.877 




Martin [7] 




4.423 


2.461 


1.763 


1.394 


1.164 


1.004 




Log [ID] 




4.706 


2.401 


1.636 


1.250 


1.015 


0.857 




Cornell [37] 




14.060 


5.681 


1.449 


3.672 


3.322 


3.088 




Buchmuller- 


-Tye M 


6.253 


3.068 


2.356 


2.032 


1.845 


1.721 




Lichtenberg 


-Wills [38] 


6.662 


3.370 


2.535 


2.139 


1.902 


1.740 
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Table 3. Mass spectra (in GeV^) of cc states. 



Meson 
State 


0.7 


Potential index 

0.8 0.9 1.0 


V 

1.1 


1.3 


Expt. 



EFG03 

US! 


ZVR95 

m 


BGE05 




i oi 


^ DQ1 


3.093 


3.095 


3.097 


^ DQQ 


O.lUl 






^ 1 nn 






2.999 


2.993 


2.987 


2.982 


2.977 


2.968 


2.980 


2.979 


3.000 


2.982 


1 3 Pn 


S 4^0 


3.488 


3.523 


3.557 


^ ^8Q 


S fi47 


O. 00\J 




S ^40 




1 3 P 1 


3.419 


3.451 


3.481 


3.508 


3.535 


3.582 


3.511 


3.510 


3.500 


3.505 


1 3 P 


3.403 


3.432 


3.459 


3.484 


3.508 


3.549 


3.415 


3.424 


3.440 


3.424 




3.434 


3.469 


3.502 


3.533 


3.562 


3.614 




3.526 


3.510 


3.516 




S 4^7 


3.519 


3.580 


3.641 


^ 700 


S 81 ^ 


o.uou 


o.uou 


^ 7^0 


^ fi79 


2 X S 


3.406 


3.459 


3.511 


3.562 


3.611 


3.707 


3.654 


3.588 


3.670 


3.630 


l 3 D 3 


3.683 


3.751 


3.817 


3.879 


3.939 


4.051 




3.815 


3.830 


3.806 


l 3 D 2 


3.694 


3.764 


3.832 


3.898 


3.960 


4.078 




3.813 


3.820 


3.800 


1 3 n 

1 


S fi77 


3.744 


3.808 


3.869 


^ Q97 


4 O^fi 


^ 770 


S 7Q8 


S 800 


^ 78^ 

J. 1 Ou 


1^2 


3.686 


3.754 


3.820 


3.883 


3.944 


4.057 




3.811 


3.820 


3.799 


9 3 Pn 


o.uuz 


3.755 


3.847 


3.939 


4 n^n 


4 9DQ 




^ Q79 


4 D9D 


1 Q79 


2 3 P X 


3.639 


3.727 


3.814 


3.900 


3.985 


4.150 




3.929 


3.990 


3.925 


2 3 P 


3.628 


3.713 


3.797 


3.880 


3.962 


4.120 




3.824 


3.940 


3.852 


9 1 P, 


O.Udl 


3.741 


3.831 


3.919 


4 DD7 


4 1 7Q 




^ Q4^ 


^ QQD 






3.673 


3.784 


3.897 


4.011 


4.125 


4.355 


4.040 


4.088 


4.180 


4.072 






3.737 


3.841 


3.945 


4 04Q 








4 1 3D 


4 OfiS 


A 3 S 1 


3.833 


3.987 


4.146 


4.309 


4.475 


4.816 


4.415 






4.406 


4^0 


3.801 


3.947 


4.097 


4.250 


4.406 


4.723 








4.384 




3.962 


4.155 


4.356 


4.564 


4.780 


5.229 










5% 


3.935 


4.120 


4.312 


4.511 


4.716 


5.140 












4.073 


4.300 


4.540 


4.792 


5.055 


5.609 










6% 


4.049 


4.269 


4.500 


4.742 


4.994 


5.522 
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Table 4. Mass spactra (in GeV) of be states. 



Meson 
State 


0.7 


Potential index 

0.8 0.9 1.0 


V 

1.1 


1.3 


AEH05 

M 


EFG03 

m 


ZVR95 

I3S] 




6.328 


6.329 


6.330 


6.331 


6.332 


6.333 


6.416 


6.332 


6.340 




6.283 


6.280 


6.278 


6.275 


6.273 


6.269 


6.380 


6.270 


6.260 


1 3 P 2 


6.700 


6.737 


6.772 


6.804 


6.835 


6.891 


6.837 


6.762 


6.760 


1 3 P 1 


6.685 


6.719 


6.752 


6.781 


6.809 


6.860 


6.772 


6.749 


6.740 


1 3 P 


6.678 


6.711 


6.741 


6.770 


6.797 


6.845 


6.693 


6.734 


6.730 


l l Pl 


6.693 


6.728 


6.762 


6.793 


6.822 


6.876 


6.775 


6.699 


6.680 


2 3 S, 


6.709 


6.770 


6.831 


6.890 


6.949 


7.063 


6.896 


6.881 


6.900 




6.684 


6.741 


6.798 


6.852 


6.906 


7.011 


6.875 


6.835 


6.850 


1 3 £> 3 


6.947 


7.016 


7.083 


7.146 


7.208 


7.322 


7.003 


7.081 


7.040 




6.952 


7.022 


7.090 


7.155 


7.218 


7.335 


7.000 


7.079 


7.030 


l 3 Di 

1 


6.944 


7.013 


7.079 


7.142 


7.202 


7.315 


6.959 


7.077 


7.020 




6.948 


7.017 


7.085 


7.148 


7.210 


7.325 


7.001 


7.022 


7.010 


2 3 P 2 


6.918 


7.011 


7.103 


7.194 


7.285 


7.462 


7.186 


7.156 


7.160 


2 3 P 


6.908 


6.998 


7.087 


7.176 


7.263 


7.435 


7.136 


7.145 


7.150 


2 3 P 


6.902 


6.991 


7.080 


7.166 


7.252 


7.421 


7.081 


7.126 


7.140 


2 l P^ 

J_ 


6.913 


7.004 


7.095 


7.185 


7.274 


7.448 


7.139 


7.091 


7.100 


3 3 5i 


6.930 


7.041 


7.154 


7.268 


7.382 


7.611 


7.215 


7.235 


7.280 


y ^ u 


6.911 


7.019 


7.127 


7.236 


7.345 


7.564 


7.198 


7.193 


7.240 


4 3 5x 


7.093 


7.247 


7.407 


7.570 


7.737 


8.079 


7.468 






4^0 


7.077 


7.228 


7.384 


7.542 


7.704 


8.035 


7.452 








7.224 


7.418 


7.620 


7.829 


8.046 


8.498 








5^0 


7.211 


7.401 


7.599 


7.804 


8.015 


8.455 








6 3 ^ 


7.337 


7.565 


7.807 


8.060 


8.324 


8.883 








6^0 


7.325 


7.550 


7.788 


8.036 


8.295 


8.841 









M(l 1 5 )=6.286 Expt. [5]. 
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Table 5. Mass spectra (in GeV) of bb states. 



Meson 
State 








Potential index 


V 




Expt. 




EFG03 

m 


ZVR95 
[39] 


0.7 


0.8 




0.9 


1.0 


1.1 


1.3 


l 3 ^ 


9.458 


9.459 


9.460 


9.460 


9.461 


9.463 


9.460 


9.460 


9.460 


1% 


9.407 


9.404 


9.402 


9.400 


9.397 


9.393 




9.400 


9.410 


l 3 Po 
x 1 2 


9.855 


9.8 


191 


9.926 


9.958 


9.989 


10.046 


9.913 


9.913 


9.860 


1 3 P 


9.841 


9.8 


!76 


9.908 


9.938 


9.967 


10.019 


9.893 


9.892 


9.870 


1 3 A 


9.835 


9.8 


168 


9.899 


9.928 


9.955 


10.006 


9.860 


9.863 


9.850 


1^1 


9.848 


9.8 


!83 


9.917 


9.948 


9.978 


10.033 




9.901 


9.880 




9.861 


9.920 


9.979 


10.037 


10.094 


10.205 


10.023 


10.023 


10.020 


2% 


9.836 


9.891 


9.946 


9.999 


10.052 


10.154 




9.993 


10.000 


1 3 £> 3 


10.102 


10 


171 


10.238 


10.302 


10.363 


10.479 




10.162 


10.150 


1 3 P 2 


10.106 


10 


.176 


10.244 


10.309 


10.371 


10.489 


10.162 


10.158 


10.150 




10.099 


10 


.168 


10.234 


10.298 


10.358 


10.473 




10.153 


10.140 


1^2 


10.103 


10 


.172 


10.239 


10.303 


10.365 


10.481 




10.158 


10.150 




10.075 


10 


.165 


10.256 


10.346 


10.435 


10.611 


10.269 


10.268 


10.280 


2 3 P 


10.066 


10 


.154 


10.242 


10.330 


10.416 


10.587 


10.255 


10.255 


10.260 


2 3 P 


10.061 


10 


.148 


10.236 


10.322 


10.407 


10.575 


10.232 


10.234 


10.240 


2 1 P 1 


10.070 


10 


.160 


10.249 


10.338 


10.425 


10.599 




10.261 


10.270 


3 3 5x 


10.083 


10 


.191 


10.301 


10.412 


10.523 


10.748 


10.355 


10.238 


10.290 




10.065 


10 


.169 


10.275 


10.381 


10.488 


10.703 




10.355 


10.370 


4 3 S 1 


10.244 


10 


.394 


10.549 


10.709 


10.871 


11.207 


10.279 






4% 


10.229 


10 


.375 


10.527 


10.682 


10.840 


11.165 








5 3 ^ 


10.373 


10 


,560 


10.757 


10.961 


11.173 


11.616 


10.865 








10.360 


10 


.544 


10.737 


10.937 


11.144 


11.576 








6 3 5i 


10.482 


10 


.703 


10.938 


11.185 


11.443 


11.991 


11.019 






6% 


10.471 


10 


.688 


10.920 


11.163 


11.416 


11.952 
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(0\QfQ\V(k,e))=f v M v e!> (8) 

where k is the meson momentum, e M and My are the polarization vector and mass 
of the vector meson. In the non relativistic quark model, the decay constant can be 
expressed through the ground state wave function at the origin ippy(0) by the Van- 
Royen-Weisskopf formula [H]. The value of the radial wave function for 0~ + , (Rp) and 
for 1 , (i?y) states would be different due to their spin dependent hyperfme interaction. 
The spin hyperfme interaction of the heavy flavour mesons are small and this can cause 
a small shift in the value of the wave function at the origin. Though, many models 
neglect this difference between (-Rp) and (RV) we account this correction by considering 

< £sd >nj 



R nJ (0) = R(0) 



1 + (SF).. 



(9) 



Mi 

Where (SF)j and < e$D >nj is the spin factor and spin interaction energy of the meson 
in the nJ state, while -R(O) and Mi correspond to the radial wave function at the zero 
separation and reduced mass of the QQ system. It can easily be seen that this expres- 
sion is consistent with the relation 

„ , . 3-Rv -\- Rp , . 

m = V A (io) 

given by [12] .Though most of the models predict the mesonic mass spectrum successfully, 
there are disagreements in the predictions of their decay constants. For example, the 
ratio was predicted to be > 1 as mp < my and their wave function at the origin 
-Rp(O) ~ -Ry(O) by most of the cases [15]. The ratio computed in the relativistic methods 
[2] predicted the ratio jr- < 1, particularly in the heavy flavour sector. The disparity of 
the predictions of these decay constants play decisive role in the precision measurements 
of the weak decay parameters as well as the spectroscopic hyperfme splitting. So, we 
reexamine the predictions of the decay constants under different potential (by the choices 
of different v) schemes employed in the present work. Incorporating a first order QCD 
correction factor, we compute, 



3|#ns(0)| 



2 



here, C 2 (a s ) is the QCD correction factor given by 



7T 



m Q - m q ln m Q _ ^ P 



(12) 

_m Q + m q 

Where 5 V = | and 5 P = 2. In the case of cc and bb systems, C 2 (a s ) becomes 1 — ^-S v ' p 
as the first term within the square bracket vanishes. Our computed values of fp and 
fv without this correction and with the correction shown in brackets up to 6S states 
are tabulated in Tables [6] -El along with available experimental results and with other 
theoretical predictions in the cases of cc, be and bb systems respectively. 
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Table 6. Pseudoscalar meson decay constant fp (MeV), Vector meson decay constant 
fv (MeV") and fp/ fv of cc states (The bracketed quantities are with QCD corrections). 



CPP;, 


IS 


2S 


3S 


4S 


5S 


6S 


fp 0.7 


365(295) 


270(218) 


230(186) 


206(167) 


189(153) 


177(143) 


0.8 


377(305) 


287(232) 


248(200) 


224(181) 


207(168) 


195(157) 


0.9 


388(314) 


302(244) 


264(214) 


241(195) 


224(181) 


211(171) 


1.0 


397(321) 


316(255) 


279(226) 


256(207) 


239(193) 


226(183) 


1.1 


404(327) 


328(265) 


292(236) 


269(218) 


253(205) 


240(194) 


1.3 


417(337) 


348(282) 


315(255) 


293(237) 


276(224) 


264(213) 


M 


335±75 












m 


292±25 












fv 0.7 


419(313) 


291(217) 


243(181) 


216(161) 


197(147) 


184(137) 


0.8 


439(327) 


314(234) 


266(198) 


238(177) 


219(163) 


204(152) 


0.9 


457(341) 


336(250) 


287(214) 


259(193) 


239(178) 


225(167) 


1.0 


473(352) 


356(265) 


308(230) 


280(209) 


260(194) 


245(182) 


1.1 


487(363) 


375(280) 


o on ( c\ a v \ 

329(245) 


300(224) 


280(208) 


264(197) 


1.3 


512(382) 


412(307) 


367(274) 


339(253) 


318(237) 


303(226) 


r 1 


416±6 


304±4 


187±8 


161±10 






1 A *1 1 

[47J 


459±28 












1 A A 1 

[44J 


/I r a 1 o O 

459±28 


364±24 


319±22 


288±18 


265±16 




[ i r 1 


551 


401 










0.7 

IV 


0.87(0.94) 


0.93(1.00) 


0.95(1.03) 


0.95(1.04) 


0.96(1.04) 


0.96(1.04) 


0.8 


0.86(0.93) 


0.91(0.99) 


0.93(1.01) 


0.94(1.02) 


0.95(1.03) 


0.96(1.03) 


0.9 


0.85(0.92) 


0.90(0.98) 


0.92(1.00) 


0.93(1.01) 


0.94(1.02) 


0.94(1.02) 


1.0 


0.84(0.91) 


0.89(0.96) 


0.91(0.98) 


0.91(0.99) 


0.92(0.99) 


0.92(1.01) 


1.1 


0.83(0.90) 


0.87(0.95) 


0.89(0.96) 


0.90(0.97) 


0.90(0.99) 


0.91(0.98) 


1.3 


0.81(0.88) 


0.84(0.92) 


0.86(0.93) 


0.86 (0.94) 


0.87(0.95) 


0.87(0.94) 





0.81±0.19 













[5] -> PDG-2006, Edwards-2001, Cvetic-2004, 

01] -> Wang-2006, [IB]->Ebert-2003. 



5. Mean Square Radii and Average quark Velocity of QQ (Qeb,c) mesons 

Apart from the decay constants, fp/v, other important properties associated with a 
mesonic state are the mean square radii (r 2 ) and the mean square velocity of the 
quark/ant iquark (f 2 ). The mean square size of the mesonic states is an important 
in the estimations of hadronic transition widths [HI 09], [50] of different QQ' systems. 
The average velocity of the quark and the antiquark within a QQ bound state are 
important for the estimation of relativistic corrections and are useful particularly in the 
NRQCD formalism as well as in the estimation of the quarkonium production rates [51J. 
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Table 7. Pseudoscalar meson decay constant fp (MeV), Vector meson decay constant 
fv (MeV") and fp/ fv of be states (The bracketed quantities are with QCD corrections). 





CPP^ 


-i o 

IS 


2S 


3S 


4S 


5S 


6S 


p 

fp 


0.7 


396(355) 


289(260) 


247(221) 


222(199) 


205(184) 


192(172) 




0.8 


412(370) 


311(279) 


270(242) 


O A C / OOA\ 

z45(zzUj 


zzo(z(J4j 


215(193) 




n n 

u.y 


4zo(oo2 J 


ool(2y ( J 


on 1 l 1 \ 
zyl(zol ) 


o^t/ o /i n\ 
zb/ (z4UJ 


zoU^zzo ) 


ZO ( (zlo j 




1 n 
l.U 


a on /ono^ 
4oy(oyo J 


oo(J(ol4 J 


01 o/oon\ 
o!2(Zo(J J 


zoy(ze>yj 


070 / O /I A \ 

z / z(z44 J 


2oy(2oo ) 




1.1 


4oU(4Uo ) 


ooo(oo(J J 


oo2(zyo J 


6i[){2(o) 


zyo^zoo J 


2ol(Zo2) 




1 o 
l.d 


ACQ l A on\ 

4oo(4zU ) 


/im txn \ 
4Ul(ooy ) 


ooy(ool ) 




oo a ( onn\ 

oo4(zyy J 


ozl(zoo J 




[15J 


4oo 














[35J 


4oU±oU 














[52J 


500 












p 

fv 


0.7 


414 (349) 


296(250) 


251(212) 


225(190) 


207(175) 


194(164) 




0.8 


/I OO /0£? /I \ 

432 (364) 


320(270) 


276(232) 


249(210) 


ooi (a ntx\ 

231(195) 


01 O/l 0/l\ 

218(184) 




n n 
U.y 


44o (o(o) 


64:2{2oo ) 


zyy^zoz J 


z /o(^zoU J 


zob^zlb J 


O/i o/on/i \ 
z4z^zU4 J 




1 n 
l.U 


40o ^oyu ) 


OOoyoVO ) 


OZZ^Z ( 1 J 


zy i ^zou j 


z i y ^zoo j 


ZOO^ZZ4 J 




1 1 
1.1 


A *7R 1 A fl 1 N 

4(0 ^4U1 J 


OOO^OZO ) 


Q/i /i foorA 
o44^zyu J 


QonV97rA 
oZU^Z / UJ 


oUZ^ZOO J 


zoy^z4o j 




1 Q 
l.o 


/1Q8 f/IOrA 
4yo ^4ZU ) 


4Zl^O00 ) 


OO ( ^oZD ) 


o04^oU ( ) 


o4o^zyo ) 


OOO^ZoZ ) 




rfrn 


quo 














[ooj 


















500 












fp 

fv 


0.7 


0.96(1.02) 


0.98(1.04) 


0.98(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.8 


0.95(1.02) 


0.97(1.03) 


0.98(1.04) 


0.98(1.05) 


0.99(1.05) 


0.99(1.05) 




0.9 


0.95(1.01) 


0.97(1.03) 


0.97(1.04) 


0.98(1.04) 


0.98(1.05) 


0.98(1.04) 




1.0 


0.95(1.01) 


0.96(1.03) 


0.97(1.03) 


0.97(1.04) 


0.97(1.04) 


0.98(1.04) 




1.1 


0.95(1.00) 


0.96(1.02) 


0.97(1.03) 


0.97(1.03) 


0.97(1.03) 


0.97(1.04) 




1.3 


0.94(1.00) 


0.95(1.01) 


0.95(1.02) 


0.96(1.02) 


0.96(1.02) 


0.96(1.02) 




ma 


0.86 













S5] 1.00 

E2] l.oo 

[15]^ Ebert-2003, [35] -»• Gerstein-1995, |52]^ Eichten-1994. 



We compute the mean square radii as 

/■OO 

<r 2 )= / r 4 |i? ni (r)| 2 rfr (13) 



and the average mean square quark/antiquark velocity for the cc and bb systems, 
according to the relation given by 



((v q ) 2 ) = ^(E-(V(r))) (14) 

Here, E is the binding energy of the system, Mi is the reduced mass of the mesonic 
system and (V(r)) is the expectation value of the potential. In the be case, the velocity 
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Table 8. Pseudoscalar meson decay constant fp (MeV), Vector meson decay constant 
fv (MeT^) and fp / fv of bb states (The bracketed quantities are with QCD corrections). 



CPP^ 


-I o 

IS 


2S 


3S 


4S 


5S 


6S 


fp 0.7 


708(606) 


492(421) 


414(355) 


370(317) 


340(291) 


318(273) 


0.8 


733(628) 


528(453) 


453(388) 


409(350) 


379(325) 


357(306) 


0.9 


756(647) 


563(482) 


490(420) 


448(384) 


419(359) 


397(340) 


1.0 


776(665) 


596(511) 


527(451) 


486(416) 


457(392) 


436(373) 


1.1 


794(680) 


627(537) 


562(481) 


523(448) 


495(424) 


A ^7 A / A /~\ /~i \ 

474(406) 


1.3 


824(706) 


686(587) 


629(539) 


594(509) 


569(488) 


549(471) 


fv 0.7 


722(584) 


497(402) 


417(337) 


372(301) 


342(276) 


320(259) 


0.8 


749(606) 


534(432) 


457(369) 


412(333) 


382(309) 


359(291) 


0.9 


773(625) 


571(462) 


495(401) 


452(365) 


422(341) 


o a a / ono\ 

399(323) 


1.0 


>-7r\ r- f n a o\ 

795(643) 


pat/ ^ oa\ 

605(489) 


roo / /|oi \ 

533(431) 


491(397) 


462(373) 


/ion/or/ 1 ^ 

439(356) 


l.i 


814(658) 


638(516) 


570(461) 


v on t a o o \ 

529(428) 


501(405) 


479(388) 


1.3 


847(685) 


Tnr\ / r f~* f~*\ 

700(566) 


641(518) 


605(489) 


578(468) 


558(451) 


[ r 


715±5 


498±5 


430±4 


336±18 


369±42 


A ,( A 1 O O 

240±28 


1 1 P I 

[16J 


839 


562 










m 


498±20 


366±27 


304±27 


259±22 


228±16 




0.7 

TV 


0.98(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.8 


0.98(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.9 


0.98(1.04) 


0.99(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


1.0 


0.98(1.03) 


0.99(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


1.1 


0.98(1.03) 


0.98(1.04) 


0.99(1.04) 


0.99(1.05) 


0.99(1.05) 


0.99(1.05) 


1.3 


0.97(1.03) 


0.98(1.04) 


0.98(1.04) 


0.98(1.04) 


0.98(1.04) 


0.98(1.04) 



[5] -> PDG-2006, [SHWang-2006, P|-> Ebert-2003. 



of 6 and c quarks are obtained as 



{{v b f) = {E- (V(r))) 2 m l (15) 

m b (m 6 + m c ) 

{(Vc f ) = (E - { V { r))) 1 ^— ) (16) 

The computed rms radii up to 6S" states of cc, 6c and bb systems are listed in Table [9] 
for the range of potential index 0.7 < v < 1.3. The estimated rms velocity (vfy 2 of the 
charm and beauty quark/antiquark using Eqn. [14] to [16] are given in Table [T0l of cc , be 
and bb systems in their IS", IP, ID and 2S to 65* states. 

6. Decay rates of quarkonia 



The spectroscopic parameters including the predicted masses and the resultant radial 
wave functions are being used here to compute the decay rates. We consider the 
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Table 9. Mean Square radii (fm) for the QQ (Q e b, c) states in various potential 
power index, v. 



CPP^ 

V 


IS 


IP 


2S 


ID 


2P 


3S 


IS 


5S 


6S 


cc 0.7 


0.50 


0.79 


1.07 


1 no 

1.03 


1.33 


1.59 


2.08 


2.55 


3.00 


0.8 


0.48 


0.74 


0.99 


0.96 


1.23 


1.46 


1.89 


2.30 


2.69 


0.9 


0.46 


0.70 


0.94 


0.90 


1.14 


1.35 


1.74 


2.10 


2.44 


1.0 


0.45 


0.67 


0.89 


0.85 


1.07 


1.26 


1.61 


1.93 


2.23 


1.1 


0.43 


0.64 


0.84 


0.81 


1.01 


1.19 


1.50 


1.79 


2.06 


1.3 


0.41 


0.59 


0.77 


0.74 


0.92 


1.07 


1.33 


1.56 


1.79 


[23] 


0.39 




0.82 






1.44 


2.36 






[54j 


0.43 




0.85 






1.18 


1.47 







be 


0.7 


0.39 


0.62 


0.84 


0.82 


1.06 


1.26 


1.65 


2.02 


2.38 




0.8 


0.38 


0.58 


0.79 


0.76 


0.97 


1.16 


1.50 


1.83 


2.14 




0.9 


0.36 


0.55 


0.74 


0.71 


0.91 


1.07 


1.38 


1.67 


1.94 




1.0 


0.35 


0.53 


0.70 


0.67 


0.85 


1.00 


1.28 


1.53 


1.78 




1.1 


0.34 


0.51 


0.67 


0.64 


0.80 


0.94 


1.19 


1.42 


1.64 




1.3 


0.33 


0.47 


0.61 


0.59 


0.73 


0.85 


1.05 


1.24 


1.42 


bb 


0.7 


0.25 


0.41 


0.55 


0.54 


0.70 


0.83 


1.10 


1.35 


1.60 




0.8 


0.24 


0.38 


0.52 


0.50 


0.65 


0.77 


1.00 


1.22 


1.43 




0.9 


0.23 


0.36 


0.49 


0.47 


0.60 


0.71 


0.92 


1.11 


1.30 




1.0 


0.22 


0.35 


0.46 


0.45 


0.57 


0.66 


0.85 


1.02 


1.19 




1.1 


0.22 


0.33 


0.44 


0.42 


0.53 


0.62 


0.79 


0.95 


1.09 




1.3 


0.21 


0.31 


0.40 


0.39 


0.48 


0.56 


0.70 


0.83 


0.94 




m 


0.19 




0.40 






0.71 


1.17 


1.85 






[53] 


0.23 
0.20 




0.51 
0.48 






0.71 
0.72 


0.88 
0.92 








m 


0.24 




0.51 






0.73 


0.93 







Vinodkumar-1999, [54J— > Gunar-1997 ,[53] -> Juan-Luis-2008, 
m\ -> Eichten-1980. 



conventional Van Royen-Weisskopf formula for the di-gamma and di-leptonic decay 
widths. Like in many other theoretical models, we also consider the contribution from 
the radiative corrections to these decays. Accordingly, the two photon decay width of 
the pseudoscalar meson is computed as 

To-+— 77 = To + Tfl (17) 
where T is the conventional Van Royen-Weisskopf formula given by 



-j^f \P^s(0)\ 2 (18) 
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Table 10. Average quark Velocity in QQ (Q e b, c) states with various potential 
power index. 



ppp 

y 


1 Q 


1 p 


9Q 
Zo 


1 Pi 


op 


oo 




Do 


Do 


cc:(v 2 c y 0.7 


0.245 


0.263 


0.268 


0.297 


0.289 


0.295 


0.320 


0.341 


0.360 


0.8 


0.265 


0.296 


0.309 


0.343 


0.339 


0.351 


0.388 


0.421 


0.450 


0.9 


0.283 


0.329 


0.350 


0.389 


0.391 


0.410 


0.462 


0.509 


0.550 


1.0 


0.300 


0.361 


0.392 


0.435 


0.445 


0.472 


0.541 


0.603 


0.660 


1.1 


0.316 


0.392 


0.434 


0.480 


0.499 


0.536 


0.625 


0.705 


0.779 


1.3 


0.345 


0.451 


0.518 


0.569 


0.610 


0.669 


0.804 


0.928 


1.043 




0.27 


0.29 


0.35 


0.34 


0.39 


0.44 


0.52 







be: 0.7 


0.396 


0.419 


0.427 


0.472 


0.459 


0.469 


0.507 


0.540 


0.570 


0.8 


0.427 


0.472 


0.492 


0.544 


0.539 


0.558 


0.616 


0.667 


0.713 


0.9 


0.457 


0.524 


0.558 


0.617 


0.622 


0.652 


0.733 


0.806 


0.871 


1.0 


0.484 


0.575 


0.624 


0.690 


0.706 


0.749 


0.859 


0.956 


1.045 


1.1 


0.509 


0.624 


0.691 


0.762 


0.793 


0.851 


0.991 


1.118 


1.234 


1.3 


0.555 


0.719 


0.825 


0.904 


0.970 


1.063 


1.276 


1.471 


1.653 



be: {vlY 


0.7 


0.030 


0.032 


0.032 


0.036 


0.035 


0.036 


0.038 


0.041 


0.043 




0.8 


0.032 


0.036 


0.037 


0.041 


0.041 


0.042 


0.047 


0.051 


0.054 




0.9 


0.035 


0.040 


0.042 


0.047 


0.047 


0.049 


0.056 


0.061 


0.066 




1.0 


0.037 


0.044 


0.047 


0.052 


0.054 


0.057 


0.065 


0.073 


0.079 




1.1 


0.039 


0.047 


0.052 


0.058 


0.060 


0.065 


0.075 


0.085 


0.094 




1.3 


0.042 


0.055 


0.063 


0.069 


0.074 


0.081 


0.097 


0.112 


0.126 





0.7 


0.077 


0.074 


0.076 


0.081 


0.079 


0.081 


0.086 


0.092 


0.096 




0.8 


0.083 


0.083 


0.087 


0.094 


0.093 


0.096 


0.105 


0.113 


0.121 




0.9 


0.088 


0.093 


0.098 


0.106 


0.107 


0.113 


0.125 


0.137 


0.148 




1.0 


0.093 


0.101 


0.110 


0.119 


0.122 


0.130 


0.147 


0.163 


0.178 




1.1 


0.097 


0.110 


0.122 


0.132 


0.137 


0.147 


0.170 


0.191 


0.210 




1.3 


0.105 


0.127 


0.145 


0.157 


0.168 


0.184 


0.220 


0.252 


0.283 




[53] 


0.094 




0.091 






0.103 


0.120 










0.080 


0.068 


0.081 


0.075 


0.085 


0.096 


0.112 







>Juan-Luis-2008, [54j^ Gunar-1997. 



and Tr is the radiative correction given by [25] 

r R = SH. (^) r (i9) 



7T \ 3 

Similarly, the leptonic decay widths of the vector mesons with radiative correction is 
computed as 

Tl — -+1+1- = ^VW + Trad (20) 
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Table 11. 0+^77 decay rates (in keV) of heavy qurkonia states (Bracketed 
quantities are with QCD corrections). 



CPP V IS 


an 


oc 
OO 


A C 


5b 


PC 1 

DO 


p 


r (r 77 ) 


-p /-p \ 

1 0(1 77) 


1 0(1 77) 


p /p \ 

1 0(1 77) 


p /p \ 

1 o(l 77) 


p /p \ 
1 o(l 77) 


cc 0.7 


5.87(3.98) 


2.83(1.92) 


1.92(1.30) 


1.47(1.00) 


1.20(0.82) 


1.02(0.69) 


O.o 


p on / a op\ 

6.29(4.26) 


3.15(2.13) 


2.17(1.47) 


1.68(1.14) 


1.38(0.93) 


1.17(0.80) 


u.y 


0.00(4.01) 


3.44(2.33) 


2.40(1.63) 


1.87(1.27) 


1 f A / 1 r\ a\ 

1.54(1.04) 


1.31(0.89) 


l.U 


0.yo(4. id) 


3.70(2.51) 


2.61(1.77) 


2.03(1.38) 


1.68(1.14) 


1.43(0.97) 


1.1 


i\20(4. yz ) 


3.93(2.67) 


2.79(1.89) 


10/1 

2.18(1.48) 


1 on/1 oo\ 

1.80(1.22) 


1.53(1.04) 


1 
l.o 


I . f4(o.24J 


a 00/0 no\ 

4.33(2.93) 


no / o no\ 

3.08(2.08) 


2.40(1.62) 


1.97(1.33) 


1.66(1.13) 


[ k 

m 


f .z±U. / 


1.3 ± 0.6* 










[00 


7 k in 
/.0 — 1U 


or a r 

3.5 — 4.5 










[06 


7.14±0.95 


AAA _i_n /I 
4.44±U.4o 










116 


5.5 


1.8 










00 U. / 


U.44(U.oo ) 


0.20(0.15) 


0.14(0.11) 


0.11(0.08) 


0.09(0.07) 


0.08(0.06) 


0.8 


0.47(0.36) 


0.23(0.18) 


0.17(0.13) 


0.13(0.10) 


0.11(0.09) 


0.10(0.07) 


0.9 


0.50(0.38) 


0.26(0.20) 


0.19(0.15) 


0.16(0.12) 


0.13(0.10) 


0.12(0.09) 


1.0 


0.53(0.40) 


0.29(0.22) 


0.22(0.17) 


0.18(0.14) 


0.16(0.12) 


0.14(0.11) 


1.1 


0.55(0.42) 


0.32(0.25) 


0.25(0.19) 


0.21(0.16) 


0.18(0.14) 


0.16(0.12) 


1.3 


0.60(0.45) 


0.38(0.29) 


0.31(0.23) 


0.26(0.20) 


0.23(0.18) 


0.21(0.16) 


W. 


0.56 


0.269 


0.208 










0.384±0.047 


0.191±0.025 










[16j 


0.35 


0.15 


0.1 









[5]-> PDG-2006 ,[55] -^Lansberg-2008, [56] -> Kim-2005, 



[16] -^Ebert-2003 ,* -^Anser-2004 [57]. 



where Tyvy is the conventional Van Royen-Weisskopf formula given by 

4o? 2 e 2 

r vw = \Rns(o)\ 2 (21) 

and the radiative correction T ra d is given by 

16 

Trad = —-^—C*s TyW, (22) 

on 

It is obvious to note that the computations of the decay rates and the radiative cor- 
rection terms described here require the right description of the meson state through 
its radial wave function at the origin R(0) and its mass M which in turn depend on 
the model parameters like a s , confinement strength and quark model masses. Gener- 
ally, due to lack of exact solutions for colour dynamics and with the uncertainties over 
the exact nature of interquark potential, R(0) and M are also been considered as free 
parameters of the theory [58] . However, we found it appropriate to employ the spectro- 
scopic parameters of the mesons such as the phenomenologically predicted meson mass 
and the corresponding wave function predicted by different models for the estimation 
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Table 12. 1 — > l + I decay rates (in keV) of heavy qurkonia states (Bracketed 
quantities are with QCD corrections). 



CPP^ 


IS 


2S 


3S 


4S 


5S 


6S 


V 


Fvw(Xu) 


rvw(r«) 




rVw(r«) 


Fvw(Xu) 




cc 0.7 


5.64(2.77) 


2.44(1.20) 


1.60(0.79) 


1.21(0.59) 


0.97(0.48) 


0.82(0.40) 


0.8 


6.19(3.04) 


2.78(1.36) 


1.85(0.91) 


1.41(0.69) 


1.14(0.56) 


0.96(0.47) 


0.9 


6.69(3.28) 


3.12(1.53) 


2.10(1.03) 


1.61(0.79) 


1.30(0.64) 


1.10(0.54) 


1.0 


7.16(3.51) 


3.45(1.69) 


2.35(1.15) 


1.80(0.88) 


1.47(0.72) 


1.24(0.61) 


1.1 


7.60(3.73) 


3.78(1.85) 


2.60(1.27) 


1.99(0.98) 


1.62(0.80) 


1.37(0.67) 


1.3 


8.38(4.11) 


4.41(2.17) 


3.07(1.51) 


2.37(1.16) 


1.92(0.94) 


1.62(0.79) 


[ r- 1 


5.55±0.14 


2.48±0.06 


0.86±0.07 


0.58±0.07 






m 


6.7(5.4) 


3.2(2.4) 










bb 0.7 


1.37(0.84) 


0.62(0.38) 


0.43(0.26) 


0.33(0.21) 


0.28(0.17) 


0.24(0.15) 


0.8 


1.47(0.91) 


0.71(0.44) 


0.51(0.31) 


0.41(0.25) 


0.34(0.21) 


0.30(0.18) 


0.9 


1.57(0.97) 


0.81(0.50) 


0.59(0.37) 


0.48(0.30) 


0.41(0.25) 


0.36(0.22) 


1.0 


1.65(1.02) 


0.90(0.56) 


0.68(0.42) 


0.56(0.34) 


0.48(0.30) 


0.43(0.26) 


1.1 


1.74(1.07) 


1.00(0.62) 


0.77(0.47) 


0.64(0.39) 


0.56(0.34) 


0.50(0.31) 


1.3 


1.88(1.16) 


1.19(0.74) 


0.95(0.59) 


0.81(0.50) 


0.71(0.44) 


0.64(0.40) 





1.34±0.018 


0.612±0.011 


0.443±0.008 


0.272±0.029 


0.31±0.071 


0.13±0.03 


m 


1.4(1.3) 


0.6(0.5) 











[5] -> PDG2006, p] -> Ebert-2003. 



of the decay properties of the mesons. 

Making use of the model parameters, the resultant radial wave functions and the mesonic 
mass we compute the 0~ + — > 7 7 and 1 — > l + l~ decay widths for each cases of the po- 
tential model employed here for the present study. The results are shown in Table (fill) 
for _+ -^77, and in Table (fl2l) for 1 — > / + l~ in comparison with the predictions 
of the contemporary potential models and with the known experimental values. The 
bracketed quantities listed in both the tables are the decay widths with the respective 
radiative corrections added to the conventional V-W formula as per EqnJTT] and Eqn. 
1201 respectively. 

7. Results and Discussion 

We have employed the coulomb plus power potential form to study the mass spectrum 
and decay properties of heavy mesons. Unlike in our earlier studies using variational 
approach [MJ EE] , here we solved the Schrodinger equation numerically using [34J . It 
helps us to study the mass spectrum of cc, be and bb mesons up to few excited states. 
Our potential parameters are fixed with respect to the centre of weight ground state 
IS" mass of the QQ (Q e b, c) systems. Our predication of the excited state of these 
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mesons for the potential index v = 0.9 to 1.3 are found to be in good agreement with the 
experimental results as well as with theoretical predictions of other models. Success of 
the present study is not only related to the numerical approach but also to the fact that 
the strength of the confinement part of the potential is made state dependent according 
to the relation — 

(n+l)3 

In Table |2l we tabulate the values of the S'-wave radial wave function at the origin, 
|-Rns(0)| 2 (in GeV^ 3 ) for the S- wave of heavy QQ systems along with other models. 
These quantities are not only essential inputs for evaluating decay constants, decay 
rates, NRQCD parameters and production cross sections for quarkonium states but also 
important for the determination of hyperfme and fine splitting of their mass spectra. We 
compared our prediction for the |i?(0)| 2 with that of Martin potential [7J, Logarithmic 
potential pi)] , Cornell potential [37] , Buchmuller-Tye potential [6] and Lichtenberg- Wills 
potential |38j . We also observe that a model independent relationship for the radial wave 
function of the be with that of cc and bb system as given by [59J 

m 2 « m 2{1 - q) \Ai\ 2g (23) 

with q = 0.35, seem to hold within 2% variation for the lower states in the potential 
range 0.5 < v < 1.5 and for higher states we find the relation hold within 5% for all 
values for v studied here. 

Our results for the decay constant of pseudoscalar meson fp, vector meson fv and their 
ratio of fp/ fv with and without the QCD corrections (given in brackets) for cc, be and 
bb mesons are listed in Tables [6] to [8] respectively from IS* to 65* states. Our results are 
compared with the available experimental values [5] and with other theoretical predi- 
cations. We could see that reduction in the fp values to about 19% in the cases of cc, 
14% in the case of bb and 10% in the case of be and reduction in the fv values to about 
25% in the cases of cc, 19% in the case of bb and 16% in the case of be are attributed 
due to the QCD correction factor. Our results for 15* state of fp for cc system is in 
good agreement with the values reported by CLEO collaboration and fv with the PDG 
average value [5]. The ratio fp/ fv without the QCD correction predicted by us lie 
between 0.87 to 0.8 in the potential range of 0.7 < v < 1.3 as against the experimental 
ratio of 0.81±0.19 [5]. The predicted values of fp for 2S to 65 states are in accordance 
with other theoretical predictions. Our results for the cc meson decay constants without 
the QCD corrections are in good agreement with the experimental data, while that for 
bb system with the QCD corrections are in accordance with the experimental results 
as well as with other model predictions. The predicted properties of the be system are 
expected to be supported by the future experimental observations. 
In Table [9j we present the mean square radii of QQ (Q e b, c) systems. Our predicted 
values are in accordance with few available predictions for cc and bb states available in 
literature. However for the be system we do not find their sizes available in literature 
for comparison. 

In Table [TUl the average quark velocity at the ground state as well as at different excited 

i _ 

states (vg) 2 of QQ (Q e b, c) systems are listed for the potential index v = 0.7 to 1.3. 
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The present results are in unit of the velocity of light. Our results for cc and bb systems 
are in accordance with the existing values reported by others [5l]up to AS states. As 
expected, the quark velocity {vf) 2 increases with higher excited states. However, it 
is also been observed that with increase in the potential index u, the quark velocity 
also increases (See Table [TOT) . It corresponds to strong binding and fast motion unlike 
the usually expected case of strong binding and slow motion. The predicted quark ve- 
locity of cc system in 6S states for the potential index 1.3 is interesting as it exceeds 
unity. Probably it may be the indication of the limit at which the cc can excite. It is 
also supported by the fact that there exist little experimental evidence for the higher 
excited states of cc systems beyond 45* level. In this potential index of 1.3 the quark 
velocity approaches the velocity of light from its 45* state onwards (0.8c) warranting the 
relativistic approaches to study this states and beyond. For the choices v < 1.3, such 
problems do not seem to be important even up to the 65" states. 

In the case of bb systems up to 6S states for all the potential index studied here suggest 
the validity of nonrelativistic treatment. The 6-quark/antiquark velocities up to 65* 
states obtained here for all the choices of the potential index 0.7 < v < 1.3 lie below 
0.3c. Thus supporting the existence of higher excited states for bb system compared 
to cc system observed experimentally. In the case of be system, we have computed the 
velocity of c-quark as well as that of the 6-quark at different excited states. The charm 
quark in be system seemed to move faster than its counter part in cc system, while 
the 6-quark in be system moves slower than that in bb system. Also, the importance of 
relativistic effects to the motion of c-quark is evident for the study of its excited states 
beyond 2S level as per the velocity predictions by the choices of power index above 
0.9. This observation in our present study also support the fact that the higher excited 
levels will be loosely bound and may not be formed to be seen experimentally. Over and 
above the predicted values of {vT) 2 would be useful in the study of the decay properties 
of QQ systems using NRQCD formalism. 

Our computed values of the di-gamma and leptonic decay widths with and without the 
radiative corrections are shown in Tables [ED and [12] respectively. Our predictions for 
cc — > 77 are in good agreement with the experimental result for the potential index 
v — 1.1 to 1.3, with out the radiative corrections. But, in the case of bb — > 77 we find 
our predictions with the radiative correction are in accordance with the values reported 
by others [56l [16] . In the case of leptonic decay widths, our predictions IV w for both 
cc and bb systems are found to be slightly over estimated in the same range of potential 
index, 1.1 < v < 1.3 and that with the radiative corrections, Fu are under estimated. 
If may be the indication of the fact that these decay of quarkonia occur not at zero 
separation of the quark and antiquark but at some finite separation. We must also look 
into the various aspects of the decay of quarkonia discussed within the NRQCD like 
formalism. We envisage such attempts for our future works. 

We further conclude here that the present study of the properties of QQ (Q e b, c) sys- 
tems based on the non-relativistic coulomb plus power potential with the power index 
ranging from 0.1 to 1.5 using numerical approach to solve the Schrodinger equation is 
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an attempt to understand the exact nature of the inter-quark potential and their pa- 
rameters that provided us the spectroscopic properties as well as the decay properties 
of the QQ system. We observe that most of the properties of the QQ systems predicted 
with the potential index in the range of 0.7 < v < 1.3 are in good agreement with the 
existing experimental results as well as with other theoretical model predictions. 
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